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Abstract 

We reexamine neutron-proton pairing as a phenomenon that should be explanable in a micro- 
scopic theory of nuclear binding energies. Empirically, there is an increased separation energy when 
both neutron and proton numbers are even or if they are both odd. The enhancement is present at 
some level in nearly all nuclei: the separation energy difference has the opposite sign in less than 
1% of the cases in which sufficient data exist. We discuss the possible origin of the effect in the 
context of density functional theory (DFT) and its extensions. Neutron-proton pairing from the 
Hartree-Fock-Bogoliubov theory does not seem promising to explain the effect. We demonstrate 
that much of the increased binding in the odd-odd system might be understood as a recoupling en- 
ergy. This suggests that the DFT should be extended by angular momentum projection to describe 
the effect. 
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FIG. 1: Neutron-proton pairing effect as seen in the neutron separation energy for N = 28 as a 
function of proton number Z. There is a consistent offset of the separation energies of odd-Z nuclei 
as compared with the average of the neighboring even-Z nuclei. 

It has been known for a long time that the nuclear binding has a mild dependence on 
the combined even-odd parities of proton and neutron numbersQ, p.mj, Q. To see 
the effect, Fig. [T] shows the neutron separation energies for nuclei with neutron number 
iV = 28 as a function of proton number Z; the separation energy is expressed in terms of the 
binding energy as S n (N, Z) = B(N, Z) — B(N — 1, Z). One sees that the neutron separation 
energies for even Z are systematically larger than the average of the separation energies for 
the neighboring odd-Z nuclei. Similar behavior is found for proton separation energies S p 
in chains of isotopes. In that case S p is greater if the number of neutrons is even than when 
the number of neutrons is odd. 

To study this behavior in more detail, we examine the separation energy differences 
S n 2 P , S P 2ni defined as the difference between the separation energy and the average for the 
two neighboring nuclei. This is 



S, 



n2p 



S n (N, Z) - {S n (N, Z + l) + S n (N, Z - 1)) /2 



(1) 
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S p2n = S P (N, Z) - {{S P {N +l,Z) + S n (N - 1, Z)) /2 

for neutrons and protons, respectively. These measures were first introduced by Jensen et 
al. 4j. With our notation, the usual measure for ordinary pairing is given by (ref. 1|, eq. 
2-92,93]) 

2A n = S n2n = S n (N, Z) - ((S n (N + 1,Z) + S n (N -1,Z)) /2 

for the neutron gap, A n , and similarly S p2p gives the proton gap. Most earlier studies of 
neutron-proton pairing used different measures for the effect. In early fits of the measured 

n n 

binding energies |3|, 15|] , the effect was parameterized as 

5 ~ mod(iV, 2) mod(Z, 2) /A. 

and attributed to an enhancement in the neutron-proton interaction. In ref. [6|, the pa- 
rameterization was changed to one have an approximate A" 2 / 3 dependence on nuclear mass 
number, 

5 = Kmod(N,2)mod(Z,2)/A 2/3 . (2) 

In ref. [3] a 9-point difference formula was proposed to describe a neutron-proton pairing 
energy. This is to be compared the 6-point difference formula we use in eq. (1). We also 
mention the shell-based mass fits of Zeldesj^, which invoke a shell-dependent term similar 
to 5. 

We find that the signs of 5 , n2p and S P 2 n are remarkably consistent across the nuclear mass 
table. Taking the data from the 2003 Audi-Wapstra mass tables[8], there are 1412 nuclei 
which have values of S n 2 P that are significant, i.e., have magnitudes larger the accumulated 
error in the experimental binding energies needed to construct the difference. Of these only 
10 nuclei had a sign for S n 2 P opposite to that seen in Fig. [U Of the 1448 measured proton 
separations S P 2 n , only 9 had the opposite sign. The nuclei with significant values of S P 2 n for 
proton separations are shown in Fig. [21 with the few opposite-sign cases shown as the black 
squares. The plot for neutron separations is very similar. There is concentration in the light 
nuclei near the N = Z line, but no obvious pattern elsewhere. 

We have also examined the dependence of the magnitude of the separation energy differ- 
ences on the mass numbers, A, of the nuclei. There is a great deal of scatter as shown in 
Fig. El but the trend is consistent with an A~ 2 I 3 dependence as in eq. (T5]). The heavy line 
shows a least squares fit to the data, \S P 2 n \ = 7.3/A 2 ^ 3 MeV. The values for \S n2p \ display a 
very similar distribution. 
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FIG. 2: Nuclei with measured proton separation energy differences S P 2 n showing the cases (black 
squares) with opposite sign from the normal. 
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FIG. 3: \S p2n \ as a function of A. The line shows the A" 2 / 3 fit eq. © with K' = 7.3MeV. 



Finally, we plot in Fig. H] all the measured nuclei, distinguishing by color those whose 
\S P 2n\ is larger or smaller than the average trend. We also find a very similar pattern for 
|>S'n2p|- There is no visible dependence on shell closures. 

We now turn to the question of what is responsible for the effect. The enhanced binding 
could arise by an increased attraction between an odd neutron and an odd proton. It 
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FIG. 4: Chart of nuclides showing the distribution of S p 2 n values that are higher (red) or lower 
(blue) than the average trend |S p2 n| = 7.3/A 2 / 3 MeV. 

could also arise by a mechanism that produced an increased binding in a nucleus with even 
numbers for both protons and neutrons. There is no way to distinguish these pictures by 
the observed systematics of the separation energy differences, because even-even and odd- 
odd binding energies contribute to the separation energy difference with the same weights. 
Both pictures are consistent with the strong similarity between the proton and neutron 
separation energy differences. Still, it is important to understand the origin of the effect if 
one is to construct accurate theories of nuclear binding based on microscopic theories such 
as the self-consistent mean field theory, also called density functional theory (DFT) 
Ordinary pairing between like particles is quite well explained by the BCS or Hartree- 
Fock-Bogoliubov extension of DFT. This suggests generalizing the HFB theory to allow 
neutron-proton pairing. Certainly, at the N = Z line neutron-proton pairing is on a same 
footing as like-particle pairing. Also, the neutron-proton interaction is stronger than the 
like-particle interaction, in free space. The special effects going on near N = Z are often 
discussed as the "Wigner energy" . It is usually parameterized in a way that does not exhibit 
a neutron-proton pairing effect away from the N = Z line and we shall consider it irrelevant 
to explain the effect. We note again that over half of the opposite-sign cases are near the 
N = Z line. There have also been limited studies of neutron-proton pairing in the HFB 
theory Typically, away from the N = Z line, condensates form in the like-particle 




sectors and prevent any pairing between neutrons and protons. We therefore doubt whether 
the effect can be explained without make some extension of the usual DFT+HFB theory. 

There is a possible mechanism that only requires a mild extension of the DFT. That is 
to exploit the higher degeneracy of states in the odd-odd nucleus to recouple the neutron 
and proton more favorably. This is easiest to understand in the situations where the mean- 
field theory approaches either the spherical shell model or the strongly deformed limit. 
Indeed, Zeldes and LiranQ] may have had this mechanism in mind in their shell-based mass 
parameterization. For the shell-model limit, consider even-even nucleus (N,Z) that has a 
spherical mean field. An added neutron goes into a spherical shell j n with an energy 6j n . 
Similarly, an added proton goes into a shell j p . When there are both added neutrons and 
protons, there is an addition neutron-proton interaction energy (j n jp\Vn P \jnjp)j depending 
on the angular momentum of the pair J. The neutron separation energies for the nuclei 
with proton numbers Z, Z + 1 , Z + 2 are, respectively, 

S n (N + 1, Z) = -e jn 

S n (N + 1,Z + 1) = ~e jn - < jnjp\V n p\jnjp >J g 

S n (N + l,Z + 2) = -e jn - < Uj 2 p ) J=0 \V np \j n (j 2 p ) J=0 > jn 

In the second equation, J g denotes the angular momentum of the odd-odd nucleus ground 
state. The last equation gives the neutron separation energy for the nucleus with two ad- 
ditional protons. Here the angular momentum coupling is determined by the three-particle 
wave function. In the spherical shell model, the two protons are coupled to angular mo- 
mentum zero in the three-particle wave function \j n {jn) J= °) ■ Standard angular momentum 
recoupling gives the neutron-proton interaction as 

Jn+Jp 

< Jn(jp) J=0 \Vnp\jn(jp) J=0 £ (2J + 1) < jjp\ \V np \ JnJp >J /(2j n + l)(2j p + 1) 

Thus, in the shell model, the energy of the odd neutron when the proton number is even 
is the (2 J + l)-weighted average over the possible neutron-proton couplings. 

This value can be estimated empirically from the spectrum of the odd-odd nuclei as the 
quantity 

Jn+Jp 

5 S = £ (2J+ l)Ej/(2j n + l)(2j p + l) (3) 

J=\jn~jp\ 
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TABLE I: Comparison of neutron-proton pair interaction energies with the recoupling model, eq. 
[3l Energies are in MeV. The quantity 5 S is defined in eq. ([3]). 
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Here Ej are measured excitation energies of the levels of the multiplet in the odd-odd 
nucleus. The quantity 6 S is thus a measure of the enhancement of the neutron separation 
energy for an odd neutron in a nucleus with an odd number of protons. 

For most odd-odd nuclei, the recoupling spectrum is difficult to determine due to the 
presence of other levels. However, near doubly magic nuclei it is often possible to make a 
spectroscopic identification fig . Some cases where we could plausibly assign the members of 
the multiplet are shown in Table [B These results are also shown in Fig. [5j The recoupling 
energy 5 S has the same order of magnitude as the separation energy differences, and also 
varies from case to case in a similar way. However, there is considerable scatter leaving room 
for other mechanism to have a role. 

If the odd-odd recoupling is dominant to produce the effect, it should be suppressed 
in deformed nuclei. The reason is that the coupling of the orbitals to the symmetry axis 
removes most of the degeneracy. Only the K quantum number remains, leaving a degeneracy 
between the two states \K P — K n \ and K p + K n in the odd-odd nucleus. To see whether 
the suppression is indeed present, we have examined the separation energies differences for 
strongly deformed even-even nuclei. We took the classification of deformed nuclei from 



ref. 



14j . which used the theoretical criterion that the static deformation of the nucleus be 
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FIG. 5: Scatter plot of 5 S , eq. ((2J), compared with S n 2 P (solid squares) and S , p 2n (circles). The 
nuclei plotted are 40 K, 48 Sc, 50 Sc, 56 Co, 60 Co, 58 Cu, 208 T1, 208 Bi, and 210 Bi. 



larger than the fluctuations about the minimum. There are 92 nuclei with measured S n 2 P 
that meet the criterion. Fitting eq. to these nuclei, we find a slightly lower value for 
K, 5.7 MeV compared to 7.3 MeV. Also a larger fraction of the deformed nuclei have very 
small values of the separation energy differences: 23% of the deformed nuclei have S n2p 
less than 3.7/A 2 / 3 MeV versus 9% for the other nuclei. The difference is not very large, 
suggesting that other mechanism beyond the recoupling effect may be needed. For example, 
configuration mixing arising from the neutron-proton interaction might depend on whether 
those nucleons are part of a pairing condensate or not. Such a mechanism would be beyond 
the usual DFT. 
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